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MULTIPLEXED HOLOGRAPHIC DATA STORAGE IN 
BACTERIORHODOPSIN 


Abstract: High density optical data storage, driven by the information revolution, 
remains at the forefront of current research areas. Much of the current research has 
focused on photorefractive materials (SBN and LiNbO,) and polymers, despite various 
problems with expense, durability, response time and retention periods. Photon echo 
techniques, though promising, are questionable due to the need for ciyogenic conditions. 
Bacteriorhodopsin (BR) films are an attractive alternative recording medium. Great 
strides have been made in refining BR, and materials with storage lifetimes as long as 
100 days have recently become available. The ability to deposit this robust 
polycrystalline material as high quality optical films suggests the use of BR as a 
recording medium for commercial optical disks. Our own recent research has 
demonstrated the suitability of BR films for real time spatial filtering and holography. We 
propose to fully investigate the feasibility of performing holographic mass data storage in 
BR. Important aspects of the problem to be investigated include various data multiplexing 
techniques (e.g. angle- amplitude- and phase- encoded multiplexing, and in particular 
shift-multiplexing), multilayer recording techniques, SLM selection and data readout 
using crossed polarizers for noise rejection. Systems evaluations of storage parameters, 
mcluding access times, memory refresh constraints, erasure, signal-to-noise ratios and bit 
error rates, will be included in our investigations. 
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INTRODUCTION 


Bacteriorhopsin(BR) provides an attractive medium for optical signal/image processing 
and optical memory applications. Attractive features of BR include its relatively low cost, ability 
to deposit it in high quality optical films on planar or convex/concave surfaces, and its low 
sensitivity to environmental factors such as temperature/humidity fluctuations. 

Our work with BR started several years ago, when M.S. student Jin Choi investigated 
holographic multiplexing techniques, using angle- and phase-encoding to multiplex multiple 2-D 
data patterns within a single BR film. 

The current NASA grant began in September of 1996, and the research has been carried 
out largely by Post-doctoral research assistant Arkady Bablumian. The objective of this research 
is to investigate BR for use in optical (holographic) memory applications. A major part of this 
investigation is to explore various holographic multiplexing techniques (angle, amplitude- shift- 
and phase-encoded- multiplexing), and model and measure various performance parameters. 
Emphasis has been placed particularly on shift-multiplexing techniques, due to their mechanical 
simplicity and effectiveness (which makes them particularly viable for commercial application). 

Our first year research results have resulted in several novel approaches, applicable not 
only to BR-based memories, but to other optical memories as well. Novel approaches include use 
of highly divergent spherical reference beams (quasi-point sources implemented using optical 
fibers), coupled with the use of differing wavelengths for hologram recording and readout (to 
improve diffraction efficiency and remedy the problem of "destructive readout"). A great deal of 
modeling and experimental work (diffraction efficiency, cross-talk, spatial resolution, etc.) has 
also been carried out on these approaches. 

Much of our first year research focused on characterizing the response of BR films, 
devising robust methods for implementing shift multiplexing, and modeling and experimental 
measurement of basic multiplexed holographic optical memory architectures. To facilitate these 
investigations, experiments employed real-time (i.e. read and write is performed simultaneously - 
and often continuously as well) holographic techniques using short-lived BR films. Although the 
real-time techniques employed in the first year's research have little or no bearing in optical 
memory systems, they facilitated alignment, testing and measurement of BR parameters, and 
enabled us to perfect the differing wavelength and point source reference beam approaches that 
will likely be employed in the final BR-optical memory demonstration system. 

In the second year, we will construct, demonstrate, model, and perform detailed 
measurements on a BR-based shift multiplexed holographic memory using longer-lived BR 
films. This will enable us to perform a myriad of systems level measurements such as access 
times, memory refresh constraints, erasure and exposure scheduling issues, signal-to-noise ratios, 
cross-talk and bit error rates. In addition, we expect to perform some investigation into the use of 
multi-layer BR media, as a means of improving the readout selectivity of multiplexed holograms 
(i.e. reducing crosstalk). 
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SUMMARY OF PROGRESS DURING YEAR ONE 

• Developed a novel shift multiplexing approach that utilizes an addressable array of point 
sources (implemented with optical fibers). This technique has also been coupled with the use of 
differing read and write wavelengths to remedy problems associated with destructive readout. 
We have also modeled this approach, and shown experimental results to agree with our theory. 

• Developed the necessary infrastructure necessary to carry out research planned for year two 

(SLM input, frame grabber output hardware, automated shutter systems for exposure scheduling, 
automated precision positioning system etc.). ° 

• Obtained and developed computer-based control software for the 2-axis precision positioning 
system that will be used for accessing the shift multiplexing demonstration system that will be 
developed in year two of the research program. 

• A visit was made to NASA Ames, February 1997, by D.J. Mehrl, J.F. Walkup and-JrChoi. A 
progress report was presented by D.J. Mehrl. 


• Two conference papers have been presented (see bibliography), and an article has been 
submitted to Applied Optics (currently in review). Copies of these papers are included in the 
attached Appendix. 


Research Plan for Year Two 

• Acquire long-lived BR films from Bend Research, and construct a demonstration setup of a 
shift-multiplexed BR-based optical memory. Construct an experimental testbed. Perform 
numerous system level measurements. 

• Continue modeling of the BR shift-multiplexed memory architecture, and compare theoretical 
predictions with experimental results. 

• If feasible, investigate use of multi-layer BR media in order to reduce crosstalk. 
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PUBLICATIONS - YEAR ONE 


Conference Papers: 

A S- Bablumian, T.F. Krile, D.J. Mehrl, and J.F. Walkup, "Recording Shift-Selective M-Type 
Volume Holograms in Bacteriorhodopsin", SPIE Proc. . Vol. 3159, SPIE Annual Technical 
Symposium, San Diego, CA, July 1997. 

A.S. Bablumian, T.F. Krile, D.J. Mehrl, and J.F. Walkup, "Grating-Type Spatial Light Modulator 
in Bacteriorhodopsin Film", presented at the 1997 Gordon Research Conference on Optical 
Signal Processing and Holography, Meriden, NH, June 29-July 4, 1997. 


Journal Papers: 

A.S. Bablumian, T.F. Krile, D.J. Mehrl, and J.F. Walkup, "M-type Thick Holograms in 
Bacteriorhodopsin Films with a High Divergence Reference Beam", Applied Optics review). 


PERSONNEL 

Faculty: 

Dr. David J. Mehrl, Principal Investigator 
Dr. Thomas F. Krile, Co-Principal Investigator 
Dr. John F. Walkup, Project Director 

Students: 

Nabeen Surestha, Lab Assistant 

Research Associate(Post-Doctoral): 

Arkady S. Bablumian 
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— Submitted to Applied Optics 
(Currently in Review) 

M-type thick holograms in Bacteriorhodopsin films with a high 

divergence reference beam 

Arkady S. Bablumian, Thomas F. Krile, David J. Mehrl, and John F. Walkup 
Department of Electrical Engineering 
Texas Tech University 
Lubbock, TX 79409 
Abstract 

The capability of using differing read and write wavelengths for reconstructing volume 
holograms recorded in a shift multiplexing geometry is analyzed and realized for M-type 
volume holograms recorded on Bacteriorhodopsin films. The intensity distribution in the 
reconstructed wave as a function of the parameters of the recording and readout beams is 
calculated. Optimal recording and retrieving geometry as well as a precise method for tuning 
the readout setup, are suggested. 

Keywords: Bacteriorhodopsin, volume holograms, shift selectivity, differing read- 
write wavelengths, holographic data storage. 

1. Introduction 

High capacity storage of optical information can be achieved by holographic recording 
throughout the volume of a thick medium [1], The strong sensitivity of thick holograms to 
the parameters of the reconstruction beam allows one to multiplex many holograms in the 
same volume and then read them out selectively. The simplest way to retrieve any desired 
hologram is to use a readout beam identical, in all respects, to the reference beam used for 



recording. At the same time, for nondestructive retrieval of stored holograms, it is necessary 
to either fix them or reconstruct them at wavelengths to which the material is insensitive. 
(The latter is of particular interest in the general area of real-time optical data processing). 
In general, changing the readout wavelength will not only change the output’s orientation 
(in accordance with the Bragg condition), but the shape of the wavefront as well. This leads 
to both a loss of information in the object beam, and to difficulties associated with synthesis 
of a readout beam with a complicated wavefront shape. The simplest case of reconstructing 
a volume hologram formed by two plane waves was considered by Kulich [2]. 
Unfortunately, from the practical point of view this approach cannot be Applied for 
analyzing several interesting hologram multiplexing geometries in situations where 
recording beams are highly divergent spherical waves with arbitrary orientations. In this 
paper we consider the results of both theoretical and experimental investigations of the 
differing read-write wavelengths reconstruction of optical information stored in a thick 
hologram with a recording geometry that is highly sensitive to a shift of the readout beam 
[1]. For an arbitrary wavelength, we calculate the optimal parameters for both the readout 
spherical beam and the appropriate intensity distribution and resolution in the reconstructed 
beam. These parameters are then compared with values derived experimentally. 

It is necessary to notice that the readout beam has to be calculated and oriented with 
respect to the hologram with high accuracy, according to the shift selectivity of a chosen 
setup [2], Since the orientation has three degrees of freedom, and recording beams have to 
be measured with the above-mentioned accuracy, it becomes clear that the results cannot be 
easily applied without some practical method of finding the optimal location of the readout 
beam. So, our next step was to develop a method for tuning the readout setup. Finally, 



experimental verification of the theory is demonstrated using M-type holograms in 
Bacteriorhodopsine (BR) films. 

2. Theory 

Let the hologram be recorded by reference (ui) and object (U2) waves as shown in 
Fig.l. In both recording arms the laser beams are expanded and then point-focused. This 
allows both the formation of a spherical reference beam and insertion of the information 
into the object arm, using a transparency with transmittance t. The scaled Fourier-transform 
of t is produced at the plane of the beam focus [(X2,y2) in Fig. 1 ]. Here we let (xi,yi) 
designate the coordinates of the spherical reference source and (X2,y2) designate the 
coordinates of the zeroth-order component of the Fourier-spectrum of transparency t. 

The hologram of transparency t can itself be regarded as the superposition of an array of 
elemental holograms, each being formed by a pair of spherical beams, whose sources are 
located at the point (xi,yi), and one of the points of the Fourier-plane (x2+Ax,y2+Ay) Fig. 1 ]. 
The hologram of the zeroth-order component of the Fourier-spectrum T [H(Ax,Ay)=H( 0 , 0 )J, 
will determine the background in the reconstructed image I, where the holograms of higher- 
order components [H(Ax,Ay)] will determine its resolution. It is evident, that each of the 
elemental holograms can be reconstructed by the beam whose parameters coincide with 
those of one of the recording beams. The existence of reconstructed beams whose 
wavelengths differ from the recording wavelength becomes clear from the following 
consideration. Each elemental hologram is recorded by point-source beams, therefore any 
of its small parts is described by the constant grating vector K. This means that for an 
arbitrary wavelength X in every point (x,y) of the elemental hologram, two rays (incident and 
diffracted) oriented at the Bragg angle to the grating vector K are determined 



unambiguously. Let us neglect the change in the grating vector K along the propagation path 
for each reconstructed ray in the hologram. In this case the lattice of rays combines to form 
only two beams which can be diffracted by this hologram and which, generally speaking, for 
an arbitrary wavelength X, are not spherical. Let us consider the hologram reconstruction by 
a spherical wave with coordinate (X3, y3) and determine the Bragg mismatch angle, 8, for its 
rays at every point of their crossing within the hologram. The assumption we have made 
above permits us to consider the diffraction of each pair of interacting rays independently 
and to calculate pointwise the diffraction efficiency of incident rays as a function of angle 5 
in accordance with coupled-wave theory [4]. Thus we can calculate the intensity'cfistribution 
in the reconstructed beam for every elemental hologram H(Ax,Ay). Their sum corresponds to 
the Fourier-spectrum of the transparency t by which we can estimate the resolution in the 
reconstructed image. 

As indicated in Fig.l, rays from the recording (X|,yO, (xa,y2) and reconstructing (X3,y3) 
sources to an arbitrary point (x,0) on the recording medium make an angle 0j where 
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which directly follows from the Bragg law written in common form [Ref.3] for readout (X 2 ) 
and recording (Xi) wavelengths. Consequently the Bragg mismatch angle S Bragg is: 


^ Bragg (*) — @ Bragg ( X ) ^i( X ) ■ 
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We choose the coordinates of the point source reconstruction beam in such a way 

that two of its rays, crossing the hologram at points (0,0) and (c,0), make Bragg angles 
with the input plane (Eq. 2.2). We then determine the mismatch angle Seragg for the rest of 
the rays of the reconstructed beam and insert it into the relationship which expresses the 
angular sensitivity of a thick hologram 



where d is the thickness of the hologram and x is the coupling constant [Ref .4]. 

This expression describes, in the exit plane of the hologram, the profile of the intensity 
of the diffracted wave which attains its maximum at points x=0 and x=c. Notice that we are 
interested in the relative diffraction efficiency of diffracted rays versus only beragg for each 
of them, so we can simplify Eq. (2.4) to obtain 

f- = V ! sinc I {*-[v I +M fcB , Wr}- < 2 - 5 > 

Here we took into account the fact that the presence of loss has very little influence on 
the angular sensitivity [Ref.4] and used Eq. (2.4) written for lossless gratings with real x- 
The constant v was measured experimentally. Figure 2 represents the relative intensity 
distribution profile (equivalent to normalized diffraction efficiency) of diffracted beams in 
the exit plane of the hologram, calculated for the elemental hologram H(0,0) [corresponding 
to the zeroth-order component of the Fourier-spectrum, T of transparency t] for 



reconstructed beams with different values of the parameter ‘c\ Values of Xi =458nm and X 2 
=633nm were used for this and subsequent plots. The X-axis is in the hologram recording 
plane and the different values of ‘c’ are obtained by changing the location of the 
reconstruction point source, (X 3 ,y 3 ). Optimized readout parameters will correspond to the 
portion of the hologram where the curves of Fig. 2 have sufficient diffraction efficiency, say 
not less then 1/e, and are as flat as possible. In our experiments, (Section 3), the hologram 
was about 3mm in extent, so c = 2mm would be optimal from Fig. 2. 

Profiles of diffracted beams with optimized parameter ‘c\ calculated for several 
recording set-ups of H(0,0) holograms (with different orientations and average spatial 
frequencies of their recording beams) are shown in Fig. 3. 

The diffraction of the beam, optimized for the reconstruction of the hologram H(0,0), 
from the holograms H(Ax,Ay), corresponding to the higher-order components of the Fourier- 
spectrum T, is shown in Fig.4. These plots represent the sum intensity distribution of beams 
diffracted by holograms H(Ax,Ay) and H(-Ax,-Ay), where we use the fact that symmetric 
components of the Fourier-spectrum contain the same information. From Fig.4, we see that 
the relative diffraction efficiency decreases with increasing values of r, corresponding to 
holograms of higher components of the Fourier-spectrum, T. Thus the reconstruction system 
is acting like a low-pass filter, and one can use this curve to define an effective spatial cutoff 
frequency, at some particular value of relative diffraction efficiency (1/e). 

Let’s consider now in more detail how the intensity profiles of beams, diffracted by 
elemental holograms (Figs. 2-4) can be utilized to determine the resolution and intensity 
distribution in the reconstructed image, t. Every point P of the recording transparency, t is 
formed by a cone of rays connecting it with points of the Fourier-spectrum (Fig.l), each of 



which we took as an elementary point-source object beam of the elemental hologram. The 
part of the transparency, t in the vicinity of the point P will be reconstructed with resolution 
h, if all rays of a cone with base radius 

r = (2.6) 

h 

fj 

have enough diffraction efficiency ( rj > ), where r=(Ax 2 +Ay 2 ) 0 5 is one of the 

coordinates of the Fourier-plane of the transparency t, and D is the distance between the 
plane of transparency t and its Fourier-transform plane, T. Taking into account the behavior 
of the plots in Fig. 4, the proper diffraction efficiency of only edge cone rays corresponding 
to the maximal-order components of the Fourier-spectrum will provide the required 
resolution. Then, it is easy to see that the intensity distribution of the reconstructed beam at 
the image plane corresponds to a profile of the beam diffracted by the H(0,0) hologram at 
its exit plane reduced in size by D 2 cos 0 av er/Di, where Di and D 2 are the distances between 
the Fourier-plane of image i and, accordingly, the hologram and the plane of the image, and 
©aver is the average incident angle of the reconstructing beam. Thus, the intensity and 
resolution in an arbitrary part of the reconstructed image T can be estimated by the profile of 
two beams diffracted by elemental holograms H(0,0) and H ( Ax , Ay ) , where 

Ax and Ay obey the condition of Eq.(2.6). 

3. Experimental Results 

To investigate the different wavelength regimes of recording and reconstructing volume 
holograms, we used a thick Bacteriorhodopsin(BR) film. This material exhibits excellent 
properties for real-time optical processing applications as well as for permanent high-density 



data storage [Ref.3]. The hologram recording in BR-films was carried out in the M-type 
regime, where the readout red light performed, at the same time, the excitation function for 
the BR molecules, transferring them into the state (M-state) ready for recording by blue 
light. M-type holograms were used as an example of nondestructive restoration, which 
simultaneously allows us to change and compare in real-time the recording and readout 
parameters. The mutant- variant BRc^N-films with a thickness of 100pm were used. The 
hologram of a standard USAF test-pattern was recorded at X=458nm and reconstructed at 
X=633nm. The setup geometry for recording and reconstruction was chosen according to the 
theory above enabling us to compare experiments results with the results derived 
theoretically (see Figs. 2 - 4). 

As was pointed out above, in the recording regime, with a high-sensitivity to shift, the 
process of finding the location of the readout beam with the given parameters becomes a 
practically difficult problem. This is particularly true for high values of the Klein parameter 
Q, which characterizes thick gratings. For optimally tuning the readout setup the following 
method was suggested and utilized. The thick recording medium was first replaced by a 
structure consisting of two separated thin films. The hologram recorded in such a structure, 
in contrast to a common thick hologram, can be reconstructed by any orientation of readout 
beams. This forms two images, propagated, in general, in different directions. Upon optimal 
orientation of the readout beam, these directions coincide and interferometric fringes of 
overlapped images coalesce into one fringe, by the maximal size of which one can judge the 
accuracy of the reconstruction source location. For such a structure, modeling the volume 
medium, a 1mm glass plate with photo resist layers on both surfaces was used. The readout 
beam was formed using a piece of optical fiber, so as to be able to manipulate the beam’s 



parameters in real-time. The setup allowed us to position the recording medium used in the 
same place and compensate for the influence of their different thicknesses on the setup 
parameters upon making the transition from the model medium to the BR-film. 

Figure 5 presents the photos of the USAF image reconstructed from the BR-film with 
different values of the Q-factor. For comparison, next to each photo are shown (Fig 5) the 
plots of intensity profiles of beams diffracted from two elemental holograms H(0,0) and 
H(Ax,Ay) recorded in the corresponding geometry. 

As one can see, the resolution (for instance 30pm) conforming to plots 2a, 2b, and 2c 
on Fig. 5 is provided all over the field of the reconstructed image for Q=100 and 300, and 
reduced in some parts for Q=750, which is supported upon analysis of the corresponding 
photos. The experimental sensitivity of the setup with Q=750 to the shift of the readout 
beam was approximately 20pm, deviating by 5%-6% from the theoretical prediction. The 
measured diffraction efficiency of the reconstructed images was approximately 4%-5%. 

4. Conclusions 

We have demonstrated a relatively simple method of calculating critical parameters for 
the process of reading & writing volume holograms using highly divergent beams of 
differing wavelengths. Experimental results were shown to agree with theoretical results. 
We considered the resolution and relative intensity distribution of the reconstructed image 
stored in the thick hologram by a high divergence reference beam which is arbitrarily 
oriented relative to the object beam. These parameters, as well as the optimal parameters of 
the readout beam, can be estimated by the intensity profiles of beams diffracted by two 
elemental holograms (recorded by a pair of spherical beams) whose simple calculation 
method has been presented above. A method for modeling a volume hologram using a 



medium consisting of two separated thin layers was also suggested. This approach allows 
one to experimentally find the precise optimal position of the readout beam. Also, the 
method makes accessible the modeling and direct visualization of the process of diffraction 
by thick holograms with arbitrary recording parameters, which could be difficult or 
impossible due to the characteristics of commonly used volume recording media (e. g. low 
life times, low diffraction efficiency, destructive readout, etc.). 

Acknowledgments 

We thank the Bend Research, Inc. for BR-films preparation. This work was supported by the 
NASA Ames Research Center, whose support is sincerely appreciated. 

References 

1. G. Barbastathis, M. Levene, and D. Psaltis, “Shift multiplexing with spherical reference 
beam”. Applied Optics 35, 2403 (1996). 

2. H.C. Kulich, “A new approach to read volume holograms at different wavelengths”. Opt. 
Com. 64, 407(1987). 

3. R.J. Collier, C. B. Burckhardt, and L. H. Lin in: Optical Holography, 128, (Acad. Press, 
New York. London, 1971) 

4. H. Kogelnik, ’’Coupled wave theory for thick hologram gratings”. Bell. Syst. Tech. J. 48, 


2909(1969). 



Fig. 1: Geometry of recording and readout rays. Here t and T denote 
the transparency and its scaled Fourier-transform respectively; 
i - reconstructed image; OF - optical fiber. 

Fig. 2: Normalized diffraction efficiency of H(0,0) versus 
hologram plane position for various ‘c’ parameters. 

Fig. 3: Normalized diffraction efficiency for cases: (xi,yi)=(7,10), 
(x 2 ,y 2 )=(6.5,10), (8,10), (5,10) - plots 1, 2, and 3,respectively. 

Fig. 4: Normalized diffraction efficiency of elemental holograms 
with various r versus hologram plane position. 

Fig. 5: Photos of transparency t from hologram H and corresponding plots 
of reconstructed beams at two elemental holograms H(r=0) - plot 1 
and H(r=lmm) - plot 2 in exit plane of hologram H for Q=100 (a). 


Q=300 (b), and Q=750 (c). 
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The purpose of our project was to increase data throughput in BR based pattern recognition systems (PRS). 
In this paper different options for building multichannel geometries for signal processing are considered. 
The following problems, for which solutions are necessary for creating such a PRS, have been investigated: 

1. Real time recording of several holograms and their selective reconstruction in a BR film with a short life 
time (less than the exposure time of one hologram) (Sec. 2) 

2. In Sec. 3 we consider the differing read-write reconstruction of optical information stored in a thick 
hologram with a recording geometry that is highly sensitive to a shift of the readout beam. For an arbitrary 
wavelength, we calculate both the optimal parameters for the readout beam and appropriate intensity 
distribution and resolution in the reconstructed beam. These parameters are then used for estimation of the 
loss of information associated with differing read-write reconstruction. 

3. In Sec. 4, multichannel correlator schemes in which the compared signals can have different wavelengths 
as well as the same wavelength have been investigated. It has also been shown that, for correlators working 
with one wavelength (for all the compared signals and the correlation signal), one input spatial light 
modulator (SLM) can be used. 


2. Hologram multiplexing in thick BR films with short lifetimes. 

To solve the problem of simultaneous recording and selective reconstruction of holograms using BR 
Films, we have used an encoding method analogous to shift multiplexing [4]. As in [4], we are using the 
strong selectivity of volume holograms to the parameters of the reference beam. However, unlike shift 
multiplexing that uses a shift of the recording medium, in our work each hologram of the input signals 
Si=Sioexp[ikj S r] is recorded by shifting the reference beam (Fig. 1 ). Each reference beam Rj=R,oexp[ikj g r] is 
formed by a piece of optical fiber OF*. The relative distance between the emitting edges of the OFi was 
chosen sufficient to suppress the cross-talk noise during the reconstruction stage. 

Coupling of laser beam with either fiber is done in node D which consists of an accusto optical deflector 
(AOD) and wide aperture microobjective Oj. Node D is adjusted so that at each moment one of the optical 
fibers OFj is excited depending on the voltage applied to the deflector. The LC SLM, driven by a computer, 
was used to input the optical signals into the setup. The readout beam was formed also using a piece of 
optical fiber so as to manipulate the beam’s parameters in real-time. The holograms have been recorded 
with wavelength X]=0.45p and reconstructed with /. 2 =0.63 p. Holograms were recorded in the M-type 
regime in which, as was mentioned above, the readout red beam performs, at the same time, the excitation 
function. The holograms have been recorded using the following method. N signals Sj are displayed on the 
SLM with frequency v. Reference beam switching is done at the same frequency. The SLM and AOD were 
synchronized so that each signal Sj from N signals S n would correspond to a single reference beam R, from 
N reference beams R„(n=l...N). Thus each pair of beams (Sj,Rj), representing the i-th hologram Hj, expose 
the recording medium with frequency v/N in the time interval 1/v. In time N/v, the recording medium will 
be exposed to all holograms H,. If this time is significantly less than the lifetime of the excited M-state of 
BR ( determining the film’s memory), then the recording conditions of the holograms are the same with 
respect to each other. To qualitatively evaluate the effect of recording N holograms on the diffraction 
efficiency of each individual hologram, we have to take into account the following: 

Let’s consider the average intensity Io of each pair (S„R,) to be the same, and pairs (S„R,) themselves as 
uncorrelated (sequential recording) and assume high spatial frequency functions (interference pattern of the 
two coherent beams Sj and Rj). In this case the total exposure E of the recording medium by all N-l 
holograms(except the i-th hologram) in each point of the recording medium is equal to - 
f V 


E=l\ 




N * 

= (A + ^ ) ~ T(N — l)l Q , where T is exposure time. 

/i=i 


(2.1) 


Thus the cumulative effect of N-l holograms on recording the i-th hologram can be considered as a 
uniform iIlumination(for large values of N) of the recording region by a beam, the intensity of which 
exceeds N-l times the average intensity of recording pairs (Rj, Sj). Figure 2a shows the diffraction 
efficiency of an M-type hologram as a function of intensity of a blue beam uniformly illuminating the 
recording region of the hologram. The plot is constructed for different values of the I 2 /Ii ratio which are 1, 
2, and 5, respectively. Figure 2b shows the similar plots of the DE versus units of the average blue intensity. 
It can be seen from Fig. 2b that by selecting a proper parameter ratio I 2 /I,, say curve 1, for small values of 



N, say, N=4 which corresponds to the (N-l)Io=3Io point, one obtains a factor of 5 decrease of diffraction 
efficiency, DE in hologram Hi compared to the single-exposure case, which will correspond to 0.5% - 0.8% 
for BR. 

In our work, simultaneous hologram recording using the above described approach has been tested for 
the case N=2. Maximum diffraction efficiency of the holograms recorded in BR was 1-1.3% instead of the 
expected value of 2%, which was evaluated using Figures 2a, b. 

Such a decrease (1.5 times ) in diffraction efficiency can be explained by the instability of the voltage 
applied to the AOD, which leads to a phase drift of R;, causing the diffraction efficiency to decrease. 


3, The influence of differing read-write wavelength reconstruction on the parameters of the 

volume holograms 

Let the hologram be recorded by reference (uj) and object (U 2 ) waves as shown in Fig. 3. In both 
recording arms the laser beams are expanded and then point-focused. This allows both the formation of a 
spherical reference beam and insertion the information into the object arm, using a transparency with 
transmittance t. The Fourier-transform of t is produced at the plane of the beam focus [ (X 2 ,y 2 ) in Fig. 3]. Let 
(x b yi) designate the coordinates of the spherical reference source and let (X 2 ,y 2 ) designate the coordinates of 
the zeroth-order component of the Fourier-spectrum of transparency t. 

The hologram of transparency t can be regarded as the superposition of an array of elemental holograms, 
each being formed by a pair of spherical beams, whose sources are located at point (x^yj), and one of the 
points of the Fourier-plane (x 2 +Ax,y 2 +Ay). The hologram of the zeroth-order component of the Fourier- 
spectrum T [H(Ax,Ay)=H(0,0)j, will determine the background in the reconstructed image I, while the 
holograms of higher-order components [H(Ax,Ay)] will determine its resolution. It is evident that each of the 
basic holograms can be reconstructed by the beam whose parameters coincide with those of one of the 
recording beams. The existence of reconstructed beams whose wavelengths differ from the recording 
wavelength becomes clear from the following consideration. Each elemental hologram is recorded by point- 
source beams, therefore any of its small parts is described by the constant grating vector K. This means that 
for an arbitrary wavelength X in every point (x,y) of the elemental hologram, two rays (incident and diffracted) 
oriented at the Bragg angle to the grating vector K are determined unambiguously. Let us neglect the change 
in the grating vector K along the propagation path for each reconstructed ray in the hologram. In this case the 
lattice of rays combines to form only two beams which can be diffracted by this hologram and which, 
generally speaking, for an arbitrary wavelength X are not spherical. Let us consider the hologram 
reconstruction by a spherical wave with coordinate (X 3 ,y 3 ) and determine the Bragg mismatch angle, 8, for its 
rays at every point of their crossing the hologram. The assumption we have made above permits us to consider 
the diffraction of each pair of interacting rays independently and to calculate pointwise the diffraction 
efficiency of incident rays as a function of angle 8 in accordance with coupled-wave theory [6]. Thus we can 
calculate the intensity distribution in the reconstructed beam for every elemental hologram H(Ax,Ay). Their 
sum corresponds to the Fourier-spectrum of the transparency t by which we can estimate the resolution in the 
reconstructed image. 

As indicated in Fig. 3, rays from the recording (xi,yj), (x 2 ,y 2 ) and reconstructing (x 3 ,y 3 ) sources to an 
arbitrary point (x,0) on the recording medium make an angle 0j where 


( 


0, = arctan 


V 


i= 1,2,3 


(3-1) 


The Bragg angle at the same point is 

6 Bragg = arcsin-i 


lx ■ 

J L V cm 

0,(x)-0 2 (-*)T| 

^ q A bill 
[^2 

2 Jj 


0l(*) + 02 (*) 


(3.2) 


which directly follows from the Bragg law written in common form [Ref. 5] for read-out X 2 and recording X\ 
wavelengths. Consequently the Bragg mismatch angle <5 Bragg ls: 

S Bragg = e Bragg( X )~ e 3( X ) ■ (3-3) 

We choose the coordinates of the point source reconstruction beam (X 3 ,y 3 ) in such a way that two of its rays, 
crossing the hologram at points (0,0) and (c,0), make Bragg angles with the input plane (Eq. 3.2). We then 



determine the mismatch angle 8 Bragg for the rest of the rays of the reconstructed beam and insert into the 
relationship which expresses the angular sensitivity of a thick hologram [Ref. 6], 


Tj = exp 


2 ad 


x exp[2i^] x 


jexp /(<f + V 2 ) 05 - exp -i(£ 2 + v 2 ) 05 ]J 


4x 


£ 2 +v 2 


(3.4) 



where d is the thickness of the hologram, % is a coupling constant. 

This expression describes, in the exit plane of the hologram, the profile of the intensity of the diffracted 
wave which attains its maximum at points x=0 and x=c. Notice that we are interested in the relative diffraction 
efficiency of diffracted rays versus only 5 Bragg for each of them, so we can simplify Eq. (3.4) to 

^ = V 2 sine 2 jtf-'fv 2 + (3.5) 

Here we took into account that the presence of loss has very little influence on the angular sensitmty [Ref. 6] 
and used Eq. (3.4) written for lossless gratings with real x* The constant v was measured experimentally. 
Figure 4 represents the relative intensity distribution profile (equivalent to normalized diffraction efficiency) 
of diffracted beams in the exit plane of the hologram, calculated for the elemental hologram H(0,0) 
[corresponding to the zeroth-order component of the Fourier-spectra T of transparency t] for reconstructed 
beams with different values of the parameter ‘c\ Values of = 0.45pm and X 2 = 0.63pm were used for this 
and subsequent plots. The X-axis is in the hologram recording plane and the different values of V are 
obtained by changing the location of the reconstruction point source, (x 3 ,y 3 ). Optimized readout parameters 
will correspond to the portion of the hologram where the curves of Fig. 4 have enough diffraction efficiency, 
say not less than 1/e, and are as flat as possible. In our experiments the hologram was about 3mm in extent, so 
c = 2mm would be optimal from Fig. 4. 

Profiles of diffracted beams with optimized parameter V calculated for several recording set-ups of 
H(0,0) holograms (with different orientations and average spatial frequencies of recording beams) are shown 
in Fig. 5. The diffraction of the beam, optimized for the reconstruction of the hologram H(0,0), from the 
holograms H(Ax,Ay), corresponding to the higher-order components of the Fourier-spectrum T, is shown in 
Fig.6. These plots represent the sum intensity distribution of beams diffracted by holograms H(Ax,Ay) and 
H(-Ax,-Ay), where we use the fact that symmetric components of the Fourier-spectrum contain the same 
information. From Fig.6, we see that the relative diffraction efficiency decreases with increasing values of r, 
corresponding to holograms of higher components of the Fourier-spectrum, T. Thus the reconstruction system 
is acting like a low-pass filter, and one can use this curve to define an effective spatial cutoff frequency, at 
some particular value of relative diffraction efficiency (1/e). 

Let’s consider now in more detail how the intensity profiles of beams, diffracted by elemental holograms 
(Fig. 4 to 6) can be utilized to determine the resolution and intensity distribution in the reconstructed image, t. 
Every point P of the recording transparency, t is formed by a cone of rays connecting it with points of the 
Fourier-spectrum (Fig.4), each of which we took as an elementary point-source object beam of the elemental 
hologram. The part of the transparency, t, in the vicinity of point P will be reconstructed with resolution h, if 
all rays of a cone with base radius 

0.611, D 

r = — (3.6) 

h 


have enough diffraction efficiency ( Tj > 


Vo 


), where r=(Ax 2 +Ay 2 ) 0 ' 5 ■ 


is one of the coordinates of the 


e 

Fourier-plane of the transparency t and D is the distance between the plane of transparency t and its 
Fourier-transform plane, T. Taking into account the behavior of the plots in Fig. 6, the proper diffraction 
efficiency of only edge cone rays corresponding to the maximal-order components of the Fourier-spectrum 
will provide the required resolution. Then, it is easy to see that the intensity distribution of the reconstructed 



beam at the image plane corresponds to a profile of the beam diffracted by the H(0,0) hologram at its exit 
plane, reduced in size by D 2 cos0 avw /Di, where Dj and D 2 are the distances between the Fourier-plane of 
image i and, respectively, the hologram and the plane of the image; 0 avcr is the average incident angle of the 
reconstructing beam. Thus, the intensity and resolution in an arbitrary part of the reconstructed image T can 
be estimated by the profile of two beams diffracted by elemental holograms H(0,0) and H(Ax, Ay) , 

where Ax and Ay obey the condition of Eq.(3.6). 

4. Multichannel correlation schemes. 

The main results of the previous two sections which can be used to advantage in building multichannel 
correlators are: 

1. The ability to simultaneously record several holograms in thick BR films that can be reconstructed 
selectively. 

2. The capability of reconstructing volume holograms using different wavelengths. The resolution of a 
reconstructed image can be sufficiently high (comparable with TV standards). The simplified method 
provided for computing the resolution and the DE of the reconstructed beam allows one to assess the 
information loss in each particular instance of recording. 

Using the above results, let’s consider one of the possible correlator setups that has several independent 
channels for processing the input signal(Fig. 7). Each of the channels is acting in a regime analogous to that 
of a CFP. Changes that have been introduced into the setup (for the reasons which will come cfear bellow) 
include: 

• In the correlator of Fig. 7 the exact Fourier-transform is changed to the scaled Fourier-transform. 

• Instead of a reference beam with a plane wavefront, we use a highly convergent one. 

Holograms H* have been recorded using optical signals Su (filter signal) formed by the input setup SLMi 
and one of the optical fibers OFj and have been used as a matched filter for the signal S 2i (master signal) 
input into the system by SLM 2 . Hologram multiplexing is done using the method described in Section 2. For 
the setup to be used as a multichannel correlator the following element adjustments should be made: 

The red beam 2 illuminating SLM 2 is formed by the lens L 2 and N pieces of fiber OF 2 and represents N 
converging spherical waves u r Parameters of each of the spherical waves are adjusted so that the diffraction 
requirements of one of the holograms, H it recorded without modulation (by SLMi t2 ) are met. In this case 
each of the spherical beams Hj will reconstruct one of the virtual images, Irrii of beam R\. The method 
applied to form a Fourier-transform allows one to change the transform’s scale by shifting SLM 2 
longitudinally along the axis of the lens L 2 , thus compensating for scale changes caused by different 
wavelengths of the signals S\ and S 2 . The SLM 2 position with respect to the axis L 2i which corresponds to 
an equal scale of the optical signals Si and S 2 , is the same for all holograms. In the direction normal to the 
optical axes of the lenses (X J 2 ) the image displayed on SLM 2 is shifted with respect to the image of SLMi, 
depending on which channel the image of SLMi has been recorded. To compensate for this shift along the 
X-axis of SLM 2 , adequate shifts are introduced along the X-axis of SLMi during the recording of the 
hologram- filters Hj. Thus the correlator of Fig. 7 functions as follows: N blue input optical signals Sn are 
applied with the frequency v to the input of SLMi, forming N hologram-fi Iters H* on the BR (details in Sec. 
2). Combined red beam 1 illuminating SLM 2 provides diffraction on any of N channels ( virtual image of 
one of the Rj ) whenever the image of SLM 2 correlates with one of the images on SLMj. It is worth noting 
that the Bragg diffraction regime of hologram-filters Hi provides a possibility of using independent 
channels, but limits the display accuracy of the optical signals. For example, for recording beams having a 
spatial frequency of (l/^ bluc ), a BR thickness of 100pm, and Di and D 2 (distance between the BR and SLMj, 
and SLM 2 , respectively) of 10 cm, the images have to be displayed with at least 0.5mm accuracy. In such a 
configuration the correlator can be only used to search for identical images. In this case one of the most 
attractive advantages of a CFP is lost, namely, the possibility of finding an image using an image fragment. 
However such a restriction can be substantially relaxed by decreasing the spatial frequencies of the beams 
which record hologram-filter Hj. Consequently, the number of channels is decreased. 

The above described correlation scheme in which signals have different wavelengths has the following 
disadvantages: 

1. Resolution is limited (Sec. 3) which in turn limits the amount of information being processed. 



2. Differences in the geometry of the optical signals Sj and S 2 caused by different wavelengths make the 
system very alignment sensitive. 

The next scheme, in which the signals Su and S 2 as well as the correlation signal Rj have the same 
wavelength(^i=0.45p), is free from the above disadvantages(Fig. 8).In addition, the setup allows one to use 
only one SLM. As in Fig. 7, the processed optical signals consist of filter signals S\\ that were recorded on 
BR as hologram-filters, and master signal S 2 that should be matched with the hologram-filters. Signals S {1 
and S 2 are alternately displayed on the SLM with the frequency 2v ( each signal is displayed with frequency 
v). N reference beams Rj are switched with frequency v and synchronized with the signals Sj using a method 
analogous to one used in Section 2. The separator K is implemented with two rotating shutters, and Sc ti 2 is 
used to separate correlation signals in the output plane(Fig. 8). Red light illuminating the BR is used for 
transferring the BR into the M-state in preparation for recording by blue light. In Fig. 9, synchronization of 
all elements composing the correlator is shown. It can be seen from the figure that the pair (Rj,Sj) which 

. . ( 2n-\ 2n^ 

records the hologram Hj exposes the BR in the time interval 


(2 n 2/t + f| 

U’ v J 


v J 


During the intervals 


the recording medium is illuminating only by signal S 2 . In the output plane the separator K 


in each moment of time is closed for beams Rj and in the interval 


(t- 


2n + 1 


it is open far the first 


diffraction orders of the signals S 2 on the hologram Hj, which represents the correlation signal between S\ t 
and S 2 multiplied by a constant factor. 


5. Conclusion 

The capability of a BR-based thick hologram for pattern recognition systems was considered. Real-time 
recording of several holograms and their selective reconstruction was demonstrated. We showed that using 
volume holograms in the Van-der-Lugt correlator with frequency plane allows one to perform multichannel 
processing of optical signals. Correlators in which the compared signals have both different and the same 
wavelength are proposed. It was shown that these setups allow one to increase data throughput significantly. 
We also suggested a simplified calculation method for estimating the loss of information due to recording 
and restoration by differing wavelengths. 
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Fig, 1 Hologram multiplexing setup: SLM is spatial light modulator; OFi are optical fibers forming reference beams 
Ri: node D consist of an acousto-optical deflector(AOD) and micro objectives Oi; BR is Bacteriorhodopsin recording 
medium 



Fig. 2: a) is relative intensity of output signal versus intensity of input signal for various values of ratio 
Ire<A>iuc (Ired/Ibiuc=5, 1, 2 for Curves 1, 2, 3, respectively), b) are the same curves given in units of 
Io(average intensity of the recording blue beams). 




X, 



Fig. 3: Geometry of recording and read-out rays. Here t and T denote the transparency and its Fourier- 
transform respectively; i - reconstructed image; OF - optical fiber 



X(mm) 


Fig. 4 Normalized diffraction efficiency of 
H(0,0) versus hologram plane position 
for various 4 c 1 parameters 


Fig. 5: Normalized diffraction efficiency for Fig 
cases: (xj.y 1 )=(7,l0). (x 2 ,y2)=(6.5.10). (8,10). 

(5,10) - plots 1,2. and 3, respectively. 





X(mm) 


Fig. 6: Normalized diffraction efficiency of elemental holograms with various r versus hologram plane 
position 



D 


Fig. 7 Multichannel correlation setup processing signals with different wavelengths: SLM, 2 are spatial light 
modulators; Lj i2 are lenses performing post-lens Fourier-transforms; 1 is a blue beam forming the object 
beam S ti for hologram Hj recording; Rj are reference beams for recording the hologram H,; 2 is a red beam 
(N spherical convergent beams Ui) which, illuminates SLM 2 , to form the master signal S 2 ; node D(acousto- 
optical modulator and N microobjectives) couples the blue laser beam with each of the OF (i ; Im Ri is one of 
the virtual images of the reference beam R ; reconstructed in red light; and DG is a diffraction grating using 
as a beam splitter. 









D 


Fig. 8 Multichannel correlation setup processing the signals with the same wavelength: SLM! is spatial light 
modulator; node D(acusto-optical modulator and N microobjectives) couples the blue laser beam with each of the OFi; 
K, consisting of revolving screens Sc 1(2 i is used for separation of the correlation signal; BR is the recording medium. 



Fig. 9 Time-dependent functions of the elements consisting of the correlator of Fig. 8: S u and R s are a filter- 
signal and a reference beam, respectively, recording the hologram, Hi (i— 1, 2,...,N); S 2 is a master-signal; Sc^ 2 
revolving shutters. 
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signal 4 ) by modulating the intensity of the input blue signal 4 . The 
blue input signal shifts the interferometric pattern into the nonlinear 
range of the characteristic curve ‘1\ thus changing the diffraction 
efficiency of the recording grating. 
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BINARY OPTICAL SIGNALS 



Fig.2: L, and L 2 are lenses, which, respectively, project (1) the input plane 
on the BR film surface and (2) the plane of the BR film onto the output plane 



Experimental setup for transformation of blue signal 
(incoherent or coherent) into the coherent red signal 







As one can see from these plots, the SLM transformation regime 
can be changed by changing Ii/I 2 . 



cr 3 c 

►— 1 CP CO 

P ™ R 

CD ^ 

_i CD O 

Sf o p 

CP rf. oo 

a < 
a « < 

P 

P CP 

CL ^ S 

P* 5* _ 

? CD O 
CP H-t 1 -^ 

HI CP 

CP y. u 

CL 

o p ttr 
p g 

r-^ A ^ ^ 

s i— i 1—1 

2 tS -1 x, 

P- HH 

o' S 'll 

g " y 

3 g* 

CP H- 1 1 
*Tl h-. ~ 
(K) J lO 

to P ^ 
w w P 5 
. t— O 

<—f H4 

H-l. ~ 

CD 

“ n 

o a 

»-b> *™t 

S- ^ 

5> 00 


HJ 

CP 

o 

o 

l-t 

CL 

i— » • 

P 

CfQ 



SJ 


to 

LO 


Fig. 3: Relative intensity of output signal versus intensity of input 
signal (input signal dan be either coherent or incoherent) for 




output signals. In our experiments these parameters were limited by the 
projection features of objective lenses LI and L2 (Fig.2) and were 
found to be 100 lines-pairs/mm for the resolution and 120 for the 
contrast. 
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Fig.4: Photos of the output plane image corresponding to the input 
transparency, T, for twd cases: (a) with the system shown in Fig. 2; 
and (b) with the BR film replaced by a mirror. 




GRAY-SCALE OPTICAL SIGNALS 



objective lenses projecting the input transparency on the BR film surface 
and the plane of the BR film onto the output plane. 
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The experimental setup for the linear regime 
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Mechanism of signal transformation in the SLM setup of Fig. 5 

Light propagation in BR is fundamentally a nonlinear process 
associated with the changing of parameters in the recording medium 
(absorption and refractive index). The blue beams used in recording the 
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Comparing this curve(Fig. 6) with the curves of Fig. 3 one can see the 
following differences between the SLM setups of Fig. 2 and Fig. 5: 



without changing the input signal’s wavelength 
it can also work in an infra-red readout mode. 
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Conclusions 



